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ABSTRACT The microbiome in fermentation has direct impacts on the quality of 
fermented foods and is of great scientific and commercial interest. Despite consider- 
able effort to explain the microbial metabolism associated with food fermentation, 
the role of the microbiome in pu-erh tea fermentation remains unknown. Here, we 
applied integrated meta-omics approaches to characterize the microbiome in two 
repeated fermentations of pu-erh tea. Metabarcoding analysis of bacterial 16S rRNA 
genes showed a decrease in the proportion of Proteobacteria and an increase in the 
abundance of Firmicutes during fermentation. Metabarcoding analysis of fungal in- 
ternal transcribed spacer (ITS) sequence demonstrated that Rasamsonia, Thermomy- 
ces, and Aspergillus were dominant at the intermediate stage, whereas Aspergillus 
was dominant at other stages in fermentation. Metaproteomics analysis assigned 
primary microbial metabolic activity to metabolism and identified microbial carbo- 
hydrate-active enzymes involved in the degradation of polysaccharides including cel- 
lulose, xylan, xyloglucan, pectin, starch, lignin, galactomannan, and chitin. Metabolo- 
mics and high-performance liquid chromatography analysis revealed that levels of 
phenolic compounds, including gallates, decreased whereas contents of gallic acid 
and ellagic acid significantly increased after fermentation (P < 0.05). The changes 
in levels of gallates and gallic acid were associated with the hydrolysis of tan- 
nase. Glycoside hydrolases, phenol 2-monooxygenase, salicylaldehyde dehydro- 
genase, salicylate 1-monooxygenase, catechol O-methyltransferase, catechol di- 
oxygenase, and quercetin 2,3-dioxygenases were hypothesized to be related to 
oxidation, conversion, or degradation of phenolic compounds. We demonstrated mi- 
crobiota in fermentation and their function in the production of enzymes related to 
the degradation of polysaccharides, and metabolism of phenolic compounds, result- 
ing in changes in metabolite contents and the quality of pu-erh tea. 


IMPORTANCE Fermented foods play important roles in diets worldwide and account 
for approximately one-third of all foods and beverages consumed. To date, tradi- 
tional fermentation has used spontaneous fermentation. The microbiome in fermen- 
tation has direct impacts on the quality and safety of fermented foods and con- 
tributes to the preservation of traditional methods. Here, we used an integrated 
meta-omics approach to study the microbiome in the fermentation of pu-erh tea, 
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which is a well-known Chinese fermented food with a special flavor and healthful 
benefits. This study advanced the knowledge of microbiota, metabolites, and en- 
zymes in the fermentation of pu-erh tea. These novel insights shed light onto the 
complex microbiome in pu-erh fermentation and highlight the power of integrated 
meta-omics approaches in understanding the microbiome in food fermentation eco- 
systems. 


KEYWORDS carbohydrate-active enzymes, food fermentation, meta-omics, 
microbiome, pu-erh tea, systemic view 


F™ ermentation is an ancient method of preserving food that depends on the repro- 
L ducible formation of multiple species of microbial communities (1). It has been 
applied to food processing for thousands of years to extend shelf life and enhance 
nutritional properties (2). It remains a major technology for essential food production 
and can be used to produce healthy foods with a wide diversity of flavors, aromas, and 
textures, including alcoholic beverages, vinegar, pickled vegetables, sausage, cheese, 
yogurt, sauces and pastes, and bread (3). Fermented foods play important roles in diets 
worldwide and account for approximately one-third of all food and beverage consumption 
(4). To date, the traditional fermentation of foods such as cider, cereal-based beverages, 
kefir, and kombucha has used spontaneous fermentation (5). Complex and dynamic 
microbial ecosystems, in which bacteria, yeast, and filamentous fungi cohabit, interact, 
and communicate, are fundamental to the quality and safety of these fermented foods. 
The microbiome also contributes to the preservation of traditional methods (1, 6). In 
addition, food fermentation is an excellent model to investigate the microbial ecology 
of natural environments, owing to its simplicity, reproducibility, and accessibility and 
the cultivability and ease of manipulation of its associated microbial communities (1, 7). 
Therefore, the microbiome in food fermentation is of great scientific and commercial 
interest. 

The field of microbiome biology (microbial community and functional analysis) has 
expanded enormously in the era of high-throughput functional genomics (8) and is 
beginning to incorporate multiple meta-omics approaches, including metagenomics, 
metatranscriptomics, metaproteomics, and metabolomics (9). Integration of these meta- 
omics approaches allows for the systematic study of microbial species, genes, metab- 
olites, and activities; it contributes to our understanding of the organization, function, 
and interactions of organisms within communities (10, 11). Meta-omics approaches 
have shed light onto complex microbial communities that were once unknown (12), for 
example, microbiota in acid mine drainage (13), gastrointestinal microbiota (14), and 
microbiota in fermentation processes (15). Therefore, integration of meta-omics ap- 
proaches is useful in the study of the microbiome in food fermentation. 

Unlike the fermentation of the most widely consumed black tea, which is actually 
the oxidation of catechins, pu-erh tea is produced by a spontaneous microbial 
fermentation. As a well-known traditional Chinese fermented food with multiple 
health benefits, such as hypolipidemic, antimutagenic, antioxidative, antitumor, anti- 
obesity, and toxicity-suppressing activities (16), it is popular in Southeast Asia and is 
becoming increasingly popular in the Western world. Microbial communities in the 
fermentation of pu-erh tea have previously been studied using culture methods and 
culture-independent approaches, including denaturing gradient gel electrophoresis, 16S 
rRNA gene clonal library, and next-generation sequencing (17). Recently, a study by Li 
et al. (18) using shotgun metagenomic and metabolomic analyses showed that signif- 
icant variations in the composition of microbiota, collective functional genes, and flavor 
compounds occurred during the fermentation process. In a previous study, we devel- 
oped an integrated metagenomics/metaproteomics approach and used the approach 
to investigate the microbial communities and enzymes of one fermented sample of 
pu-erh tea (19). However, to our knowledge, there have been no reports conducted 
using an integrated meta-omics approach to study the microbiome in pu-erh tea 
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FIG 1 Outline of this integrated meta-omics study. 


Sensory evaluation 


fermentation. In addition, little is known about the enzymes involved in the metabolism 
of phenolic compounds, which are essential for tea processing. 

To conduct a systematic review of the fermentation mechanism of pu-erh tea, the 
microbiome in fermentation was studied using an integrated meta-omics approach. 
Bacterial 16S rRNA gene and fungal internal transcribed spacer 1 (ITS1) amplicon sequenc- 
ing was used to characterize microbial succession during fermentation; the microbial 
activity was studied using a liquid chromatography-tandem mass spectrometry (LC- 
MS/MS)-based metaproteomics analysis. Metabolic succession was detected by high- 
performance liquid chromatography (HPLC) and an ultrahigh-performance liquid chro- 
matography (UPLC)-MS/MS-based metabolomics analysis. Further, we surveyed 
microbial enzymes involved in the degradation of polysaccharides and the metabolism 
of phenolic compounds (Fig. 1). This work advanced the production of pu-erh tea and 
highlights the power of integrated meta-omics approaches in understanding the 
microbiome in food fermentation ecosystems. 


RESULTS AND DISCUSSION 
Fermentation of pu-erh tea and sample collection. Two repeats of spontaneous 
fermentation of pu-erh tea (FB and FM) were developed based on the traditional 
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TABLE 1 Samples of tea leaves collected in fermentation of pu-erh tea and the designed analyses? 


Sample identifier by repeat and sample 


Collection( date, Repeat FB Repeat FM 
Sample day/mo (2014) First sample Second sample First sample Second sample 
Raw material 10/10 FBO-1 FBO-2 FMO-1 FMO-2 
Repiling 
First 18/10 FB1-1 FB1-2 FM1-1 FM1-2 
Second 25/10 FB2-1 FB2-2 FM2-1 FM2-2 
Third 30/10 FB3-1 FB3-2 FM3-1 FM3-2 
Fourth 10/11 FB4-1 FB4-2 FM4-1 FM4-2 
Fifth 20/11 FB5-1 FB5-2 FM5-1 FM5-2 
Sixth 26/11 FB6-1 FB6-2 FM6-1 FM6-2 
Seventh 29/11 FB7-1 FB7-2 FM7-1 FM7-2 
Final fermented tea leaves 1/12 FB8-1 FB8-2 FM8-1 FM8-2 


“All samples were analyzed by sensory evaluation, HPLC, the spectrophotometric method, metabarcoding analysis, and metaproteomics analysis. Additionally, raw 


material, the fourth repiling, and the final fermented tea leaves were subjected to metabolomics analysis. 


method. In total, 36 samples were collected and analyzed, as outlined in Fig. 1 and 
Table 1. Sensory evaluation revealed that the infusion of raw tea material was slightly 
astringent in taste and yellow, but after fermentation, the tea infusion became mellow 
in taste and reddish-brown. This sensory quality change was similar to general fermen- 
tation observations of pu-erh tea. It is believed that changes in sensory qualities are 
closely related to changes in chemical metabolites and caused by microbial activity. 
Therefore, the composition and activity of microbes in the fermentation process were 
further investigated by integrated meta-omics approaches. 

Microbiota in the fermentation of pu-erh tea. A total of 877,307 bacterial 16S 
rRNA and 1,740,425 fungal ITS sequences were obtained from 36 samples of tea leaves, 
and most of the microbial diversity was captured by this analysis (see Fig. S1A and B in 
the supplemental material). As the fermentation process progresses, bacterial richness 
and diversity increase, showing a peak in B4 or M6 and decreasing slightly in the last 
stages of fermentation (Fig. 2A). Fungal richness and diversity decreased in the initial 
stage of fermentation, and increased slightly during the middle and final stages of 
fermentation (Fig. 2B). 

Principal-component analysis (PCA) based on the relative abundance (RA) of oper- 
ational taxonomic units (OTUs) revealed that both bacterial and fungal communities in 
the raw material (BO and MO) differed from those of other fermented stages (Fig. 2C and 
D); bacterial communities in the final stages of the fermentation of tea leaves (B7, B8, 
M7, and M8) were similar to each other and yet they differed from those of other stages. 
However, fungal communities in the final fermentation steps of tea leaves (B8 and M8) 
were similar to those in other fermented samples (B1, B2, B3, M1, M2, M3, and M4). 

Bacterial OTUs were classified into 31 phyla (Data Set S1, sheet 1). At the phylum 
level, there was a decrease in the proportion of Proteobacteria and an increase in the 
abundance of Firmicutes. Additionally, the RA of Actinobacteria increased to 15.37% 
(B4), 24.90% (B8), and 57.51% (M7) (Fig. S1C). More than 80% of bacterial OTUs were 
assigned to 201 families (Data Set S1, sheet 2). As the fermentation process progressed, 
proportions of Enterobacteriaceae decreased in FB, whereas they increased in FM. The 
abundance of Bacillaceae increased at the final stages in both FB and FM with RAs of 
44.92% (B8) and 89.05% (M8). Additionally, RAs of Lactobacillaceae were greater than 
70% in B5 and B6, and the RA of Pseudomonadaceae was 46.56% in M1 (Fig. 2E and 
Fig. S1E). Similarly to our previous findings (20), the RA of Enterobacteriaceae decreased and 
the proportions of Bacillaceae increased during the fermentation process of pu-erh tea. 

Fungal OTUs were classified into four phyla and unclassified fungi. The dominant 
phylum was Ascomycota, which accounted for 84.76% to 99.94% of the total fungi 
(Fig. S1D and Data Set $1, sheet 3). More than 87% of fungal reads were classified into 
166 genera; Aspergillus and Debaryomyces were dominant genera in BO with RAs of 
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FIG 2 Microbiota in the fermentation of pu-erh tea. (A and B) Bacterial (A) and fungal (B) Chaol and abundance-based coverage estimates (ACE) and 
Shannon indices. (C and D) PCA of relative abundance (RA) of bacterial (C) and fungal (D) OTUs in each sample. (E and F) Circos plots of changes of major 
bacterial families (E) and fungal genera (F) (RAs, >1% of the total sequences) in fermentation. “Others” are composed of the family or genus showing 
a percentage of reads of <1.0% of the total reads in each sample. 
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55.69% and 28.90%, respectively, whereas the RA of Aspergillus was 80.03% in MO. In 
combination with Aspergillus (RAs of 10.05 to 63.78%), Rasamsonia (19.58 to 58.68%) 
and Thermomyces (3.61 to 18.7%) were the dominant fungi in the samples of B4, B5, B6, 
M5, M6, and M7, whereas Aspergillus was the dominant fungus in other fermented 
samples with RAs greater than 65% (Fig. 2F, Fig. S1F, and Data Set S1, sheet 4). The 
dominance of Aspergillus and the presence of Rasamsonia and Thermomyces in the 
fermentation of pu-erh tea have been previously reported (19, 21). 

To obtain a measure of microbial association, three OTU cooccurrence networks 
were constructed. In the bacterial cooccurrence network, most of the OTUs correspond- 
ing to Bacillaceae, Comamonadaceae, Pseudoalteromonadaceae, Pseudomonadaceae, 
Phyllobacteriaceae, and Vibrionaceae cooccurred with others in fermentation. OTUs 
assigned to unclassified groups in Rickettsiales and Enterobacteriaceae were negatively 
correlated with other bacteria (Fig. 3A). In the fungal network, most of the OTUs 
assigned to Thermomyces and Thermoascus have negative correlations with other 
genera, whereas the OTUs corresponding to Aspergillus, Rasamsonia, Penicillium, and 
Debaryomyces cooccurred with each other (Fig. 3B). The network of OTUs of 16S RNA 
genes and ITS sequences showed that bacteria and fungi were mutually exclusive in 
fermentation. For example, members of the Bacillaceae, Pseudoalteromonadaceae, Co- 
mamonadaceae, and Enterobacteriaceae showed negative correlations with fungi in 
Aspergillus, Rasamsonia, Penicillium, Debaryomyces, and Saccharomycetes (Fig. 3C). 

Overview of metaproteomics results. After protein extraction, LC-MS/MS analyses, 
and a database search, 68 to 1,582 microbial proteins in each repeated analysis and 
4,623 and 6,505 unique proteins in FB and FM, respectively, were identified and further 
annotated, with proteins identified as malate dehydrogenase, superoxide dismutase, 
and catalase, among others (Data Set S2, sheets 1 and 2). In the Gene Ontology (GO) 
annotation, the majority of identified proteins categorized as molecular functions were 
primarily catalytic activity and binding; those categorized as biological process were 
cellular process and metabolic process; the proteins categorized as cellular components 
were cell parts and protein-containing complex. The majority of enzyme classes were 
oxidoreductases, transferases, and hydrolases (Fig. S2A and Data Set S2, sheets 3 and 
4). The major categories identified by the Cluster of Orthologous Group analysis were 
energy production and conversion, translation, ribosomal structure and biogenesis, 
posttranslational modification, protein turnover and chaperones, amino acid transport 
and metabolism, and carbohydrate transport and metabolism (Fig. S2B and Data Set S2, 
sheet 5). A total of 116 KEGG pathways were annotated, with the most common 
pathways identified as glycolysis/gluconeogenesis (ko00010), ribosome (ko03008), ox- 
idative phosphorylation (ko00190), and citrate cycle (ko00020) (Fig. S2C and Data Set 
S2, sheet 6). These KEGG pathways grouped into cellular processes, environmental 
information processing, genetic information processing, and metabolism. KEGG path- 
ways in metabolism were further associated with classes of amino acid metabolism, 
biosynthesis of other secondary metabolites, and carbohydrate metabolism (Fig. S2C). 
Overall, the majority of identified proteins were assigned to cellular process and 
metabolic process in the GO analysis and enriched in pathways belonging to metab- 
olism or genetic information processing. These data support the findings from micro- 
bial growth and reproduction. 

Metabolic succession in fermentation. The concentrations of 16 characteristic 
components of tea were measured by HPLC or spectrophotometric methods. We 
observed three change trends among the results (Fig. 4A; Table 2): (i) the levels of tea 
polyphenols (TPs), free amino acids (FAA), (—)-epigallocatechin (EGC), (+)-catechin (Q), 
1,4,6-tri-O-galloyl-B-p-glucose (GG), (—)-epicatechin 3-O-gallate (ECG), and (—)- 
epigallocatechin 3-O-gallate (EGCG) decreased with the development of fermentation; 
(ii) the contents of water extractions (WE), kaempferol, quercetin, myricetin, and 
(—)-epicatechin (EC) increased at the initial stage and then decreased significantly after 
fermentation (P < 0.05); and (iii) the levels of soluble sugar (SS), gallic acid, and ellagic 
acid increased significantly after fermentation (P < 0.05). Additionally, the caffeine 
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FIG 3. Cooccurrence networks of OTUs. (A) Bacterium-bacterium cooccurrence network. (B) Fungus-fungus cooccurrence network. (C) 
Bacterium-fungus cooccurrence network. Nodes correspond to OTUs, and connecting edges indicate negative (green) or positive (red) 


correlations between them. 


content increased significantly (P < 0.05) in FB but showed no significant change in FM 
(P > 0.05). 

A total of 11,423 m/z were detected in the metabolomics analysis (Data Set S3, 
sheets 1 and 2), and the PCA with 67.0% variation showed that the metabolites 
identified in samples from raw material (BO-1, BO-2, MO-1, and MO-2), the middle stage 
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of fermentation (B4-1, B4-2, M4-1, and M4-2), and the final stage (B8-1, B8-2, M8-1, and 
M8-2) of fermentation were distinct (Fig. 4B). Two hundred ninety-eight metabolites for 
which relative levels changed significantly in comparisons of B4 and BO, B8 and B4, B8 
and BO, M4 and Mo, M8 and M4, and M8 and MO (variable importance in projection 
[VIP] > 1.0, P<0.05, and fold change [FC] > 1.5 or FC < 0.66) were identified, for 
example, rutin, catechin 3’-O-gallate, theogallin, and (+)-catechin. The majority of 
these metabolites belong to classes of flavonoids (117 metabolites), glycerophospho- 
lipids (56 metabolites), fatty acyls (18 metabolites), carboxylic acids and derivatives (17 
metabolites), and organo-oxygen compounds (15 metabolites) (Data Set S3, sheets 2 
and 3). Among the metabolites, 44 were annotated into 61 KEGG pathways, with highly 
represented pathways of flavonoid biosynthesis (ko00941), flavone and flavonol bio- 
synthesis (ko00944), and biosynthesis of phenylpropanoids (ko01061), which included 
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eight, five, and five metabolites, respectively. Twenty-eight KEGG pathway terms were 
annotated by both metaproteomics and metabolomics, including citrate cycle 
(ko00020) and glycerophospholipid metabolism (ko00564) (Data Set S3, sheet 4). 

Relative levels of 96 to 134 or 35 to 64 metabolites were decreased (VIP > 1.0, 
P<0.05, and FC < 0.66) or increased significantly (VIP > 1.0, P< 0.05, and FC > 1.5) 
in comparison with B4/BO, B8/B4, B8/BO, M4/M0, M8/M4, and M8/MO, respectively 
(Fig. 4C). After fermentation, the relative levels of 124 and 125 metabolites were 
decreased significantly in comparisons of B8 and BO and of M8 and MO, respectively. 
Among them, relative peak areas of 64 and 57 metabolites, respectively, decreased by 
more than 10-fold, including EGCG, theaflavin digallate, luteoliflavan, and t-theanine, 
whereas the relative peak areas of 55 and 52 metabolites significantly increased after 
fermentation in comparison of B8 and BO and of M8 and MO, respectively, including 
margrapine A, kukoamine A, Thr-Trp-OH, Phe-Lys-OH, uridine, ellagic acid, and gallic 
acid (Data Set S3, sheet 5). HPLC determination verified the increasing contents of gallic 
acid, which increased 11.63 and 5.94 times in FB and FM, respectively. 

CAZyme analysis. Carbohydrate-active enzymes (CAZymes) are involved in the 
assembly and breakdown of complex carbohydrates, including oligosaccharides or 
polysaccharides as well as glycoconjugates to nucleic acids, proteins, lipids, polyphe- 
nols, and other natural compounds (22). They are responsible for the synthesis (through 
glycosyltransferases [GTs]), degradation (glycoside hydrolases [GHs], polysaccharide 
lyases [PLs], carbohydrate esterases [CEs], and auxiliary activities [AAs]) and recognition 
(carbohydrate binding module [CBM)]) of all the carbohydrates on Earth (23). Through 
dbCAN annotation, 284 and 471 proteins in FB and FM, respectively, and a total of 558 
unique proteins were annotated to 131 items in the CAZy database; the identified 
proteins were distributed into all six families of GHs, GTs, AAs, CEs, PLs, and CBMs. 
Highly represented family items were AA3 (28 proteins), AA7 (20 proteins), AA2 (15 
proteins), CBM (12 proteins), CE10 (17 proteins), GH13 (21 proteins), GH28 (18 proteins), 
GH3 (19 proteins), GH6 (23 proteins), GH72 (16 proteins), GT30 (17 proteins), GT4 (40 
proteins), and PL1_4 (6 proteins) (Fig. S3A and Data Set S2, sheet 7). 

In these CAZymes, 82 and 137 enzymes belonged to 31 and 39 families or subfam- 
ilies of CAZymes, respectively, which were hypothesized to degrade plant and fungal 
polysaccharides, including cellulose, xylan, xyloglucan, pectin, starch, lignin, galacto- 
mannan, and chitin (Fig. S3B and Data Set S2, sheet 8). For example, glucanase (GH5, 
-12, -16, -7, and -55; AA9 and AA11), cellobiohydrolase (GH7), glucosidase (GH1, -3, -17, 
and -55), and glucose-methanol-choline (GMC) oxidoreductase (AA3) may be involved 
in the biodegradation of cellulose. Peroxidase, aryl-alcohol dehydrogenase, GMC oxi- 
doreductase, choline oxidase, and alcohol oxidase, belonging to AA2 or AA3, are involved 
in the degradation of lignin. A series of enzymes, including alpha-L-arabinofuranosidase 
(GH51, -54, and -62), alpha-t-rhamnosidase (GH78), alpha-n-arabinofuranosidase (GH51 and 
-54), beta-galactosidase (GH35 and -42), galactanase (GH53), endopolygalacturonase 
(GH28), exo-a-L-1,5-arabinanase (GH93), pectate lyase (PL1_7 and -3_2), pectin lyase 
(PL1_4), polygalacturonase (GH28), rhamnogalacturan acetylesterase (CE12), rham- 
nogalacturonase (GH28), and 1,4-beta-xylosidase (GH3), are able to degrade the pectin 
backbone and the side chains of the hairy regions of pectin. Galactosidase (GH35, -42, 
-27, and -36), arabinofuranosidase (GH51, -54, and -62), xylosidase (GH3), xylanase 
(GH10, -11, and -43), and alpha-t-fucosidase (GH29) are able to degrade the xylan and 
xyloglucan. Alpha-glucosidase (GH13 and -31) and alpha-amylase (GH13) can biode- 
grade starch. Galactosidase (GH27, -35, -36, and -42), glucosidase (GH1, -3, -17, and -55), 
beta-mannanase (GH5), and beta-mannosidase (GH2) degrade galactomannan. Chiti- 
nase (GH18), glucanase (GH5_9, -16, and -17), glucanosyltransferase (GH16, -17, and 
-72), glucan-1,3-beta-glucosidase (GH17), and beta-hexosaminidase (GH20) are in- 
volved in the degradation of chitin. Therefore, the microbiota involved in fermentation 
produce CAZymes related to the degradation of plant or fungal biomass polymers and 
use the resulting monomers and oligomers for their growth and reproduction (Fig. $3B 
and Data Set S2, sheet 8). Degradation of plant polysaccharides can destroy the cells of 
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tea leaves. This is supported by the observation of the softening of tea leaves during 
the fermentation of pu-erh tea. Similarly, Wang et al. showed that the surfaces of tea 
leaves were covered by microbiota and the cell structures were largely disrupted after 
fermentation (24). 

In the metabolomics analysis, 76 glycosides were identified, with the majority being 
glucosides (31 metabolites), rhamnosides (11 metabolites), glucuronides (nine metab- 
olites), and glucopyranosides (seven metabolites). After fermentation, relative levels of 
36 and 33 glycosides were decreased (VIP > 1.0, P< 0.05, and FC < 0.66), whereas six 
and five of them showed increased levels (VIP > 1.0, P< 0.05, and FC > 1.5) in B8/BO 
and M8/MO, respectively (Data Set S3, sheet 6). GHs (EC 3.2.1.—) (25) and GTs (26) (EC 
2.4.x.y) catalyze the hydrolysis and formation of the glycosidic bond; thus, they may be 
involved in the degradation or synthesis of glycosides. For example, glucosidase (e.g., 
AOA2G7G438 and AOA1F2K440), beta-glucuronidase (AOA100IUA9), and 1,3-beta- 
glucanosyltransferase (AOAOF4YKO2) were suggested to be involved in the metabolism 
of glucosides and glucopyranosides. Galactopyranoside and galactoside were hydrolyzed 
by beta-galactosidases (AOAOF4Z133, AOAOL1J957, AOA1J6WUZ7, and AOA1L9RJP3). Alpha- 
l-rhamnosidase (AOAOR1I3B8) was suggested to hydrolyze rhamnoside (Table S1). 

Metabolism of phenolic compounds in the fermentation of pu-erh tea. Phenolic 
compounds possess one or more aromatic rings with one or more hydroxyl groups and 
generally are categorized as phenolic acids, flavonoids, coumarins, and tannins (27). 
Phenolic compounds, primarily catechins, are the characteristic chemical component in 
teas and provide a number of health benefits, such as reducing the incidence of 
coronary heart disease, diabetes, and cancer (28). Moreover, the oxidation of catechins 
and the production of oxidation reaction products in fermentation are crucial to the 
quality of black tea (29). Therefore, understanding the metabolism of phenolic com- 
pounds is essential for the investigation of the process and quality control of tea. 

In this metabolomics study, 144 phenolic compounds were identified, including 
catechin 3-O-gallate, gallocatechin, and quercetin (Data Set S3, sheet 7). After fermen- 
tation, the relative levels of 73 metabolites decreased, which was in accordance with 
the decreasing levels of polyphenols and catechins, including EGC, EC, EGCG, C, and 
ECG, as shown by the spectrophotometric or HPLC analyses. Thus, phenolic compounds 
were actively metabolized, and the majority of levels of phenolic compounds decreased 
during fermentation. The decreasing content of phenolic compounds is responsible for 
the transformation of taste from astringent to mellow. 

Seventy-two identified phenolic compounds grouped as phenolic glycosides, including 
luteone 7-glucoside, and quercetin 3,7-diglucuronide. The majority of the phenolic glyco- 
sides showed decreasing concentrations, except cerarvensin, catechin 3'-glucoside, 
terbutaline-1-glucuronide, and pelargonidin 3-arabinoside (Data Set S3, sheet 6). As 
described above, phenolic glycosides were suggested to be degraded by GHs, such as 
glucosidase, glucuronidase, galactosidase, and rhamnosidase (Table S1). 

Twenty-one identified phenolic compounds were gallates, including EGCG, epiaf- 
zelechin 3-O-gallate, and theaflavin digallate, which showed decreasing contents dur- 
ing fermentation (Data Set S3, sheet 8). For example, the content of EGCG decreased 
from 47.3 + 2.42 mg/g (BO) and 46.02 + 2.81 mg/g (MO) to 0.08 + 0.15 mg/g (B8) and 
was finally undetectable (M8). In contrast, the level of gallic acid increased more than 
11 times (Table 2). Tannase (EC 3.1.1.20) catalyzes the breakdown of ester and depside 
linkages in hydrolyzable tannins, such as tannic acid, producing gallic acid and glucose 
(30) (Fig. 5A). In the metaproteomics analysis, six tannases (AOA229XTVO, AOA124BYX3, 
AOA1L9N9KO, AOA254TWP2, AOA254U742, and AOA117DX77) and three tannases/feru- 
loyl esterases (G3YCQ1, AOA117E377, and AOA1L9N5Y7) were identified. Therefore, 
tannases were hypothesized to hydrolyze gallates and release gallic acid, which results 
in the changes of levels of these compounds in fermentation. Additionally, three 
chlorogenic acid esterases (EC 3.1.1.42) (AOA100I6W6, AOA146FON8, and AOA100IE01) 
were identified, which were suggested to hydrolyze epigallocatechin 3-O-caffeate, 
epigallocatechin 3-O-p-coumarate, and epigallocatechin 3-O-cinnamate (Fig. 5B). 
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Catechol is a phenol containing a benzene ring with adjacent hydroxyl groups. 
It can be considered the basic structural unit of phenolic compounds. According to 
the GO and KEGG annotations, we identified 19 enzymes (Table S2) involved in 
oxidation, conversion, and degradation of catechol and constructed a network as 
follows: catechol O-methyltransferase (AOAOF4YT61) catalyzes the formation of 
catechol from o-methoxyphenol, phenol 2-monooxygenase (AOAOK8LGT6) oxidizes 
phenols to form catechol, salicylaldehyde dehydrogenase (AOAOF4YRG3) oxidizes sali- 
cylaldehyde to produce salicylate, salicylate 1-monooxygenase (AOAOF4YJI4) oxidizes 
salicylate to form catechol, and catechol 2,3-dioxygenases (M2ZX00) and catechol 
1,2-dioxygenase (AOAOK8LRX4, AOA117DW12, AOA117DZB1, and AOA117DW92) de- 
grade catechol to form cis,cis-muconic acid or 2-hydroxymuconate semialdehyde 
(Fig. 5C). We hypothesized that aromatic rings or catechol in phenolic compounds was 
oxidized, modified, or degraded by these enzymes in fermentation. Additionally, quer- 
cetin 2,3-dioxygenase (AOA10011V8), which catalyzes the decyclization of quercetin, 
was identified and suggested to degrade quercetin and quercetin glycosides (Fig. 5D). 

We hypothesized a metabolic pathway of tea phenolic compounds resulting in the 
decrease of relative levels of most phenolic compounds and an increase in the content 
of several compounds including gallic acid, ellagic acid, quercetin, and myricetin in 
pu-erh fermentation as follows: (i) phenolic glycosides were hydrolyzed or synthe- 
sized by GHs and GTs; (ii) gallates were hydrolyzed by tannase and produced gallic 
acid; (iii) phenolic compounds were oxidized, modified, or degraded by catechol 
O-methyltransferase, phenol 2-monooxygenase, salicylaldehyde dehydrogenase, salic- 
ylate 1-monooxygenase, catechol 2,3-dioxygenases, catechol 1,2-dioxygenase, and 
quercetin 2,3-dioxygenase. To our knowledge, this is the first report on the enzymes 
involved in the metabolism of phenolic compounds in pu-erh tea fermentation, which 
are characteristic compounds in tea and are responsible for the taste and health 
benefits. Additionally, phenolic compounds are ubiquitously distributed phytochemi- 
cals found in most plant tissues and are important for the quality of plant-based foods 
(31); therefore, the findings in this article may provide interesting insight into other 
plant-based fermented foods and beverages. 

Development of the FFMP. Fermented foods are important societal traditions and 
are crucial regional products in terms of the economy, as well as being rich in 
microbiological resources awaiting exploration (32). Understanding the microbiome 
within fermentation ecosystems is essential for maintaining traditional and artisanal 
practices in the context of urbanization, designing starter cultures, directing sensory 
quality, and improving the safety of the consumable products (6). Previously, Parente 
(33, 34) developed FoodMicrobionet, which provides a wealth of information on the 
structure of food biomes. We suggest developing the Food Fermentation Microbiome 
Project (FFMP) to study the microbiome within the food fermentation ecosystem using 
the powerful integration of meta-omics approaches. This work provides an example of 
a study of microbiomes in a fermented food ecosystem using integrated metabarcod- 
ing, metaproteomics, metabolomics, and HPLC approaches. 

Conclusion. Microbiomes in two fermentations of pu-erh tea were systematically 
examined via the integration of metabarcoding, metaproteomics, and metabolomics 
analyses. We identified the microbial succession and association, microbial activity, and 
changes in the metabolites during the fermentation of pu-erh tea. We found that micro- 
biota produced CAZymes to degrade plant or fungal polysaccharides for their growth 
and reproduction, as well as enzymes involved in hydrolysis, oxidization, modification, 
or degradation of phenolic compounds (Fig. 6). This study advanced our understanding 
of the fermentation mechanism of pu-erh tea related to the microbial and metabolic 
succession, as well as the microbial functions during the fermentation of pu-erh tea. 


MATERIALS AND METHODS 

Fermentation of pu-erh tea and sample collection. Two traditional fermentation processes of 
pu-erh tea were developed by the Yunnan D Tea Co., Ltd., Yunnan, China, between 10 October and 1 
December 2014 (Table 1). Sun-dried green tea leaves purchased from Bajiazia Village (24°13’16.10’N, 
98°24'0.18"E) and Mingzhi Mountain (24°09'33.74"N, 98°36'51.58”E), Mangshi City, Yunnan Province, 
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Enzymes 


Salicylate 1-monooxygenase * 
Salicylaldehyde dehydrogenase « 
Phenol 2-monooxygenase .« 
Tannase + 
Glycosyl transferase * 

Catechol O- astinlansterase 
Chlorogenate hydrolase « 
Beta-glucuronidase « 
Quercetin 2,3-dioxygenase « 
Catechol 2,3-dioxygenase + 
Catechol 1,2-dioxygenase « 
1,4-beta-glucanosyltransferase « 
Oligo-1,6-glucosidase « 
Glucosidase « 
Alpha-glucosidase « 
Alpha-L-rhamnosidase « 
Endo-beta-1,4-galactanase , 
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FIG 6 Overview of microbial enzymes involved in degradation of polysaccharides and metabolism of phenolic 


compounds. 
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China, were used as the raw material. The fermentation of pu-erh tea was developed according to the 
traditional method of spontaneous fermentation; raw materials, water, utensils, and the environment 
were not sterilized, and no starter was used. Based on the temperature of the tea piles and the 
experience of the manufacturer, the tea masses were broken down, washed with water to a moisture 
content of approximately 40%, mixed, and restacked in piles for 3 to 10 days. Samples of tea leaves were 
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collected from the tea piles at five time points before each round of breaking up, mixing, and repiling. 
Samples of tea leaves were divided into two parts. One part was air dried and subjected to sensory 
evaluation, according to the protocol described by GB/T 23776-2009 (35), and analysis of the chemical 
compounds, while the second part was stored at —80°C. In total, 36 samples were collected and analyzed 
as outlined in Fig. 1 and Table 1. Detailed approaches are described in Text S1 in the supplemental 
material. 

Metabarcoding of bacterial 16S rRNA gene and fungal ITS sequence. To analyze the taxonomic 
composition of the bacterial and fungal communities, the universal primer pairs 515F (5'-GTGCCAGCM 
GCCGCGGTAA-3’) and 907R (5’-CCGTCAATTCMTTTRAGTTTI-3’) and ITS1F (5’'-CTTGGTCATTTAGAGGAAG 
TAA-3’) and ITS1R (5’-GCTGCGTTCTTCATCGATGC-3’), which incorporate Illumina adapters and barcode 
sequences, were used to amplify the V4-V5 hypervariable region of bacterial 16S rRNA genes, as well as 
the ITS1 of fungal 18S rRNA genes using a two-step amplification procedure. DNA extraction, PCR, and 
Illumina MiSeq sequencing (2- by 150-bp reads) were performed by TinyGene Technology Co., Ltd. 
(Shanghai, China). Each sample was extracted for two replicates, and each extraction was analyzed twice. 
Analysis of OTU cooccurrence networks was developed using the CoNet application (36) on Cytoscape 
3.7.1 (37). The detailed approaches are described in Text S1 in the supplemental material. 

Metaproteomics experiments. The microbial proteins in tea leaves were extracted by Tris-HCl- 
phenol and methanol precipitation, measured using the Bradford method with bovine serum albumin as 
a standard, and validated with sodium dodecyl sulfate-polyacrylamide electrophoresis (as described in 
our previous report [19]). For each sample of tea leaves, three independent extractions were carried out. 
A total of 200 yg of protein was digested with trypsin according to the filter-aided sample preparation 
protocol (38). Liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis of each replicate 
of peptide extracts was performed using an Easy-nLC1000 coupled to a QExactivePlus mass spectrometer 
(Thermo Fisher Scientific, Bremen, Germany). Raw data were processed using Thermo Proteome Discov- 
erer software version 1.4 (Thermo Fisher Scientific, Bremen, Germany) with the default settings. The 
MS/MS data were queried against the UniProt database (http://www.uniprot.org/) with the following 
search parameters: carbamidomethylation of cysteine as the fixed modification, oxidation of methionine 
and deamidation of glutamine and asparagine as variable modifications, a maximum of two missed 
cleavages, a precursor ion mass tolerance of 10 ppm, and an MS/MS tolerance of 0.05 Da. Decoy database 
searches were performed with a false-discovery rate (FDR) cutoff of 1%. GO annotations for the identified 
proteins were assigned according to those reported in the UniProt database. COG annotations of 
identified proteins were computed using eggNOG-Mapper based on eggNOG 4.5 orthology data (39, 40). 
The CAZymes annotation was developed by dbCAN (41). The KEGG pathway was annotated using the 
KEGG (Kyoto Encyclopedia of Genes and Genomes) Automated Annotation Server (KAAS) using the 
bidirectional best hit BLAST method (https://www.genome,jp/tools/kaas/) (42). Detailed approaches are 
described in Text S1 in the supplemental material. 

Metabolomics experiments. Metabolites were extracted from the raw material (BO-1, BO-2, MO-1, 
and MO-2), the fourth repiling (B4-1, B4-2, M4-1, and M4-2), and the final fermented tea leaves (B8-1, B8-2, 
M8-1, and M8-2) and assessed using a UPLC-quadrupole time of flight (Q-TOF) MS-based metabolomics 
approach. The tea powders (50 mg) mixed with 20 jl of L-2Cl-Phe (0.03 mg/ml) as internal standard were 
extracted with 1 ml of 70% methanol for 30 min in an ultrasonic bath. The extraction was kept at —20°C 
for 20 min. Then, the samples were centrifuged at 14,000 rpm at 4°C for 10 min, and 200 pl of superna- 
tant was filtered using 0.2-14m polytetrafluoroethylene filters and subjected to UPLC-Q-TOF-MS analysis 
at Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China). Triplicate preparations and analyses were 
performed for each sample. Detailed approaches are described in Text S1 in the supplemental material. 

Analysis of chemical compounds in tea leaves by HPLC and spectrophotometry. The contents 
of water extractions (WEs), tea polyphenols (TPs), free amino acids (FAA), and soluble sugar (SS) in tea 
leaves were analyzed using the spectrophotometric method described by Wang et al. (24). The amount 
of gallic acid, caffeine, hydrolyzable tannins (1,4,6-tri-O-galloyl-B-p-glucose [GG]), and catechins, includ- 
ing (+)-catechin (C), (—)-epicatechin (EC), (—)-epigallocatechin (EGC), (—)-epicatechin 3-O-gallate (ECG), 
and (—)-epigallocatechin 3-O-gallate (EGCG), in the tea leaves was determined using HPLC with an 
Agilent 1200 series HPLC system (Agilent Technologies, Santa Clara, CA, USA), as described in our 
previous study (19). Three replicates of each sample were extracted, and each extraction was detected 
twice. Data were analyzed and statistical analyses were performed in SPSS 19.0. Significant differences 
between two groups were noted by different letters (P < 0.05). Detailed approaches are described in Text 
S1 in the supplemental material. 

Data availability. The sequencing data of bacterial 16S rRNA genes and the fungal ITS1 of 18S rRNA 
gene are available at the Sequence Read Archive under project code SRP139059. The mass spectrometry 
proteomics data have been deposited in the ProteomeXchange Consortium (http://proteomecentral 
.proteomexchange.org) via the iProX partner repository (43) with the identifier PXD012223. 
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